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. The production of a Higgs boson in association with a pair of tt quarks will play a very 

(^-^ ' important role at both hadron and lepton colliders. We review the status of theoretical 

f^*) . predictions and their relevance to Higgs boson studies, with particular emphasis on the recently 

£f) ■ calculated NLO QCD corrections to the inclusive cross section for pp, pp — > tth. We conclude 
by briefly discussing the case of exclusive bbh production and the potential of this process in 

*"^h , revealing signals of new physics beyond the Standard Model. 

a: 

Q" 1 ' 1 Introduction 

Present and future colliders will play a crucial role in exploring the nature of the electroweak 
symmetry breaking and its relation to the origin of fermion masses. The discovery of a Higgs 
boson is therefore among the most important goals of both the Tevatron and the Large Hadron 
Collider (LHC). Had such a particle to be discovered, a high energy Linear Collider (LC) will 
be able to identify it unambiguously. 

The present lower bounds on the Higgs boson mass from direct searches at LEP2 are > 
114.4 GeV (at 95% CL)Hfor the Standard Model (SM) Higgs boson (h), and M h o >91.0 GeV 
and M A o > 91.9 GeV (at 95% CL, 0.5 < tan/3 < 2.4 excluded) d for the light scalar (h°) and 
pseudoscalar (A ) Higgs bosons of the minimal supersymmetric standard model (MSSM). At the 
same time, global SM fits to electroweak precision data imply Mh < 211 GeV (at 95% CL)^, 
while the MSSM requires the existence of a scalar Higgs boson lighter than about 130 GeV. 
The possibility of a Higgs boson discovery in the mass range near 115-130 GeV thus seems 
increasingly likely. 

In this context the associated production of a Higgs boson with a pair of ti quarks is kine- 
matically accessible, has a very distinctive signature, and can give the only handle on a direct 
measurement of the top quark Yukawa coupling, perhaps the most crucial coupling in exploring 
the origin of fermion masses. 

Observing — > tth at the Tevatron (y/s = 2 TeV) will require very high luminosity 1 ^ and will 
probably be at the edge of the machine capabilities. On the other hand, if < 130 Ge V, pp 
tth is an important discovery channel for a SM-like Higgs boson at the LHC (^/s = 14 TeV)^ 
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Given the statistics expected at the LHC, pp — > tih, with h — ► bb, t + t~, W + W~, 77 will also 
be instrumental to the determination of the couplings of a discovered Higgs boson E3 8 | 9 | 10 | 
Several analyses show that precisions of the order of 10-15% on the measurement of the top 
quark Yukawa coupling can be obtained with integrated luminosities of 100 fb _1 per detector. 
Morever, the combined measurements of pp — ► tih with h — > bb and h — ► t + t~ could provide the 
only model independent determination of the ratio of the bottom quark to the r lepton Yukawa 
couplings®. 

At a LC, the top quark Yukawa coupling can be measured in a model independent way via 
e + e~ — ► ti/i. The inclusive cross section for e + e~ — ► ii/i (and bbh) has been calculated including 
the first order of QCD corrections, both in the SM and in the MSSM™ and the theoretical 
uncertainty is reduced in this case to less than 10%. However, the precision of the me asu rement 
is severely limited by the machine center of mass energy. Dedicated studies show thatHH, at the 
optimal center of mass energy of -y/s^SOO GeV, integrated luminosities of 1000 fb _1 will allow 
to determine the top quark Yukawa coupling at the 5% level, for Mh = 120 GeV. However, at 
a center of mass energy of y / s = 500 GeV the e + e~ — > tih event rate is tiny and, for the same 
range of Higgs masses and integrated lumino sity, a LC will initially measure the top Yukawa 
coupling with precisions of at best 20%1 -^ l ffi Given this intrinsic limitation, the role played by 
the LHC and, in this context, by the associated production of a Higgs boson with a pair of ti 
quarks becomes even more important. 

In view of its phenomenological relevance, a lot of effort has been recently invested in im- 
proving the stability of the theoretical predictions for the hadronic inclusive total cross section 
for pp, pp — > tih. Since the tree level or Leading Order (LO) cross section is affected by a very 
large renormalization and factorization scale dependence, the first order of QCD corrections 
have been calculated and the Next-to-Leadi ng (NLO ) cross section, for a SM Higgs boson, has 
been obtained independently by two groups ^^^ l The NLO cross section has a drastically 
reduced renormalization and factorization scale dependence, of the order of 15% as opposed 
to the initial 100% uncertainty of the LO cross section, and leads to increased confidence in 
predictions based on these results. 

The calculation of the NLO corrections to the hadronic process pp, pp — ► tih presents chal- 
lenging technical difficulties, ranging from virtual pentagon diagrams with several massive in- 
ternal and external particles to real gluon and quark emission in the presence of infrared sin- 
gularities. A general overview of the techniques developed and employed in our calculation are 
presented in Section [21 and the corresponding results are illustrated in Section 03 We conclude 
with a brief outlook in Section 0J 

2 QCD corrections to tih production at the Tevatron and the LHC 

The inclusive total cross section for pp — ► tih at 0(a^) can be written as: 

Vnlo(p ( P -> tih) = V - —z ( dxxdx2 Iff {x\, fi)J^ (p \x 2 , n)o% to (xx, x 2 , fi) + 



ij 



(1) 

where J-f^ are the NLO parton distribution functions (PDFs) for parton % in a (anti)proton, 
defined at a generic factorization scale fif = fj,, and a^ LO is the 0(a^) parton- level total cross 
section for incoming partons i and j, made of the channels qq,gg — ► tih and (q,q)g — ► tih{q,q), 
and renormalized at an arbitrary scale [i T which we also take to be [i r = \x. We note that the effect 
of varying the renormalization and factorization scales independently has been investigated and 
found to be negligible. The partonic center of mass energy squared, s, is given in terms of the 
hadronic center of mass energy squared, s H , by s = xix 2 s H . At the Tevatron center of mass 
energy the cross section is entirely dominated by the qq initial state and the results presented 



in Section are obtained by including only qq — ► tih at the parton level. At the LHC center of 
mass energy the cross section is dominated by the gg initial state, but the other contributions 
cannot be neglected and are included in our calculation. 

We write the NLO parton- level total cross section x 2 , A 4 ) as: 

a^ LO (x 1 ,x 2 ,fi) = a^ (x 1 ,x 2 ,^) +da^ LO (x 1 ,x 2 ,^) , (2) 

where cr l l (x\, x 2 , fi) is the 0(a 2 s ) Born cross section, and 5a % ^ LO {x\,x 2 , n) consists of the 0(a s ) 
corrections to the Born cross sections for gg, qq — ► tih and of the tree level (q, q)g — ► tih(q, q) 
processes, including the effects of mass factorization. Sa l ^ LO {xi, x 2 , n) can be written as the sum 
of two terms: 

SaH LO {xi, x 2 , n) = J d(PSz)Y)A mT Sj - tih)\ 2 + J d{PS A )Y)A real {ij - tih + /)| 2 

= *Srt (»1 1 »2 , + ^riaJ ( x l » x 2, /*) , (3) 

where \A vir t(ij -» tt/i)| 2 and |Aeaz(u -> £i/i + /)| 2 (for ij = qq,gg and / = or ij = qg,qg and 
l = c L->q) are respectively the O(a^) terms of the squared matrix elements for the ij — ► tiTi and 
ij — > tih+l processes, and indicates that they have been averaged over the initial state degrees 
of freedom and summed over the final state ones. Moreover, d(PS3) and d(PS^) in Eq. Q denote 
the integration over the corresponding three and four-particle phase spaces respectively. The 
first term in Eq. © represents the contribution of the virtual one gluon corrections to qq — ► tih 
and gg — > tih, while the second one is due to the real one gluon and real one quark/antiquark 
emission, i.e. qq,gg — ► tih + g and qg(qg) — ► it/i + q(q)- 

The 0(ag ) virtu al and real corrections to qq — > tih and gg — > tih have been discussed in 
detail in RefeJHUSI and we will highlight in the following only the most challenging tasks. 



2.1 Virtual correction 

The calculation of the 0(a s ) virtual corrections to qq,gg — > tih proceeds by reducing each 
virtual diagram to a linear combination of tensor and scalar integrals, which may contain both 
ultraviolet (U V) an d infrared (IR) divergences. Tensor integrals are further reduced in terms of 
scal ar i ntegrals The finite scalar integrals are evaluated by using the method described in 
Ref.QSland cross checked with the FF packagel^. The scalar integrals that exhibit UV and/or 
IR divergences are calculated analytically. Both the UV and IR divergences are extracted by 
using dimensional regularization in d=4— 2e dimensions. The UV diverge nces ar e then removed 
by introducing a suitable set of counterterms, as described in detail in R e f s LL4lE] The remaining 
IR divergences are cancelled by the analogous singularities in the soft and collinear part of the 
real gluon emission cross section. 

The most difficult integrals arise from the IR-divergent pentagon diagrams with several 
massive particles. The pentagon scalar and tensor Feynman integrals originating from these 
diagrams present either analytical (scalar) or numerical (tensor) challenges. We have calculated 
the pentagon scalar integrals as linear combinations of scalar box integrals using the method 
of Ref.ELTand cross checked them using the techniques of Ref.^^l Pentagon tensor integrals 
can give rise to numerical instabilities due to the dependence on inverse powers of the Gram 
determinant (GD), GD=det(pj-pj) for pi and pj external momenta, which vanishes at the bound- 
aries of phase space when two momenta become degenerate. These are spurious divergences, 
which cause serious numerical difficulties. To overcome this problem we have calculated and 
cross checked the pentagon tensor integrals in two ways: numerically, by isolating the numerical 
instabilities and extrapolating from the numerically safe to the numerically unsafe region using 
various techniques; and analytically, by reducing them to a numerically stable form. 



2.2 Real correction 



In computing the 0{a s ) real corrections to qq,gg — > tih and (q,qg — ► ii/i + (q,q) it is crucial 
to isolate the IR divergent regions of phase space and extract the corresponding singularities 
analytical ly. We achieve this by using the phase space slicing (PSS) method, in both the doubkJ221 
and singleksSlMI cutoff approaches. In both approaches the IR region of the tih + g phase space 
where the emitted gluon cannot be resolved is defined as the region where the gluon kinematic 
invariants: 

s ig = 2pi ■ pg = 2EiE g (l - fa cos 9 ig ) (4) 
become small. Here pi is the momentum of an external (anti)quark or gluon (with energy 



Ei), f3i = Jl — mf/Ef, p g is the momentum of the radiated final state gluon ((anti)quark) (with 
energy E g ), and 9{ g is the angle between pi and p g . In the IR region the cross section is calculated 
analytically and the resulting IR divergences, both soft and collinear, are cancelled, after mass 
factorization, against the corresponding divergences from the 0(a s ) virtual corrections. 
The single cutoff PSS technique defines the IR region as that where 

S'ig ^ Smin ■> (5) 

for an arbitrarily small cutoff s m i n . The two cut-off PSS method introduces two arbitrary 
parameters, 5 S and 5 C , to separately define the IR soft and IR collinear regions according to: 

„ „ 8 8 y/s 

Eg < — - — soft region , 

(1 — cos 9i g ) < 5 C collinear region . (6) 

In both methods, the real contribution to the NLO cross section is computed analytically 
below the cutoffs and numerically above the cutoffs, and the final result is independent of these 
arbitrary parameters. With this respect, it is crucial to study the behavior of a NLO in a region 
where the cutoff (s) are small enough to justify the analytical calculations of the IR divergent 
contributions to the real cross section, but not so small as to cause numerical instabilities. 



3 Results for tth production at hadron colliders 

The impact of NLO QCD corrections on the tree level cross section is summarized in Figs.^|I]for 
both the Tevatron and the LHC. Results for a LO are obtained using the 1-loop evolution of a s (fi) 
and CTEQ4L parton distribution functions'^, while results for a NLO are obtained using the 2- 
loop evolution of a s (fi) and CTEQ4M parton distribution functions, with a^. LO {Mz) = 0.116. 

Figs, ^and |21 illustrate the renormalization/factorization scale dependence of a LO and a NLO 
at the Tevatron and the LHC. In both cases the NLO cross section shows a drastic reduction 
of the scale dependence with respect to the lowest order prediction. Figs. Eland 0] complement 
this information by illustrating the dependence of the LO and NLO cross sections on the Higgs 
boson mass at both the Tevatron and the LHC. 

The overall uncertainty on the theoretical prediction, including the errors coming from parton 
distribution functions and the top quark mass (which we take to be mt = 17 A GeV), is reduced 
to only 15-20%, as opposed to the 100-200% uncertainty of the LO cross section. Including 
NLO QCD corrections decreases (Tevatron) or increases (LHC) the LO cross section for a broad 
range of commonly used renormalization and factorization scales (obtained e.g. by varying [i 
by a factor of two around fx = fj,o), and over the entire Higgs boson mass range considered in 
our study. This can be summarized by defining a X-factor, K = cr NLO /a LO , which is however 
affected by the same strong scale dependence as the LO cross section, as well as by the choice 
of PDFs. When using CTEQ4 PDFs the K-factor corresponding to Figs. PHI is around 0.7—0.95 
at the Tevatron and 1.2 — 1.4 at the LHC, for most choices of scales and Higgs boson mass. 
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Figure 1: Dependence of <Jlo,nlo{pp — * tth) on 
the renormalization/factorization scale fi, at y / s7T = 
2 TeV, for M h = 120 GeV. 
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functions of Mh, at ^[s~h = 1 TeV, for mi. 
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Figure 3: Dependence of <j lo ,nlo{pp — * tth) on 
the renormalization/factorization scale fi, at ^J~Sh = 
14 TeV, for M h = 120 GeV. 
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4 Outlook 



The techniques developed to calculate the NLO cross section for pp,pp — ► tth can now be applied 
to the study of the associated production of bbh. The inclusive cross section for bbh production 
receives contributions from bb — > h, bg — > bh, and gg — > bbh, in order of decreasing cross section 
(qq — * bbh is negligible at both the Tevatron and the LHC). On the other hand, the exclusive 
cross section, corresponding to the experimental situation when both final state b quarks are 
tagged, receives contributions from gg — > bbh only and can be directly calculated, including NLO 
corrections, from the corresponding results for gg — * tth. In spite of the smaller cross section, 
the exclusive process is experimentally very interesting since it corresponds to a well defined 
maesurement, where final state b jets are isolated via cuts on the transverse momentum of the 
b and b quarks. The cross section for gg — > bbh is negligible in the SM and the detection of a 
Higgs boson in this channel would unambiguously signal the presence of new physics responsible 
for an anomalously large bottom quark Yukawa coupling, like the MSSM. This could actually 
be a unique opportunity within the kinematical reach of the Tevatron. 
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